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or  to  speed  calculational  time  by  code  vectorization.  In 
addition,  graphic  techniques  were  developed  to  improve  the 
display  of  the  results,  which  is  important  both  to  detect 
coding  errors  and  to  elucidate  physical  mechanisms.  In  this 
report-,  we  describe  the  codes  that  were  improved  and  present 
examples  of  the  graphical  displays  now  available  for  analysis 
of  the  calculational  results. 
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I. 


INTRODUCTION 


The  primary  objectives  of  this  research  effort  were  to  provide 
numerical  and  computational  support  in  terms  of  code  conversions,  code 
improvement,  and  the  development  of  graphic  techniques  and  codes  as  well  as 
provide  analytical  support  in  a  number  of  radiation/matter  interaction 
investigations.  Each  of  these  items  are  individually  discussed  below 
starting  with  the  various  codes  and  their  functions.  Also  a  short 
description  of  each  code  is  presented  for  the  purposes  of  documentation. 

The  Plasma  Radiation  Branch  of  NRL  is  involved  in  a  numoer  of  highly 
diverse  activities  ranging  from  studies  in  basic  atomic  physics  and 
collision  processes  to  radiation  hydrodynamics  of  imploding  targets  and 
nuclear  weapons  environment  simulations.  In  order  to  describe  the  physics 
of  such  a  variety  of  disciplines  a  number  of  computer  simulations  and  codes 
must  be  constructed  and  developed.  Most  of  these  codes  are  generated 
internally  in  the  Branch  although  some  of  them  are  transferred  from  other 
research  Institutes  and  different  computer  facilities.  It  is  inevitable 
that  some  of  these  codes  are  unvectorizable  due  to  their  logical  structure 
and  numerical  techniques  employed. 

When  a  computer  program  is  vectorized,  the  use  of  tl*a  ASC,  which  is  a 
vector  machine,  is  recommended.  In  this  case  the  ASC  runs  are  at  least  a 
factor  of  3  more  cost  effective  than  identical  runs  on  similar  non  vector 
machines.  The  VAX,  which  is  a  non  vector  machine,  is  superior  to  the  ASC 
only  when  the  codes  are  completely  unvectorized  but  again  this  superiority 
is  impaired  by  the  complete  absence  of  batch  program  facilities  and  fast 
printer  availability. 

II.  ATOMIC  PHYSICS  CODES 

A  set  of  atomic  physics  codes  from  the  NRL  Plasma  Radiation  Branch  was 
converted  from  ASC  to  VAX.  The  distorted  wave  scattering  code  "DIRECT", 
the  Dipole  Transition  Potential  TRANP  using  known  wave  functions  of  inner 
electrons,  and  the  program  "SPINFL"  which  calculates  collision  strength,  n  , 
for  any  transition  between  levels  A,  B  of  an  atom  with  2  electrons  outside 
closed  shells  have  been  improved  and  converted  to  VAX. 
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charge  of  the  nucleus, 

number  of  equivalent  core  electrons  given  by  the  wave  function 
(£  qt  =  total  number  of  core  electrons), 

adjustable  constant, 


=  dipole  polarizability  of  the  core, 
r1  =  parameter  for  the  polarization  potential. 

The  first  three  terms  on  the  right-hand  side  of  (2)  represent  the 
potential  due  to  the  nuclear  charge  and  the  screening  of  core  electrons 
(spherically  averaged).  The  fourth  term  is  the  "exchange-correlation 
potential"  similar  to  the  definition  given  by  Slater  (but  not  identical  to 
it),  and  the  last  term  is  the  polarization  potential. 

This  form  of  V(r)  is  convenient  if  of  core  electrons  are  known. 
Parameters  and  r^  for  the  polarization  potential  are  known  only  rarely 
and  may  be  omitted. 
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(b)  Analytical  potential 


V(r)  =  -  -  Z  +  (Z  -  Z)  exp  (-cBr) 
**  n 


where 

Zn  =  charge  of  the  nucleus, 

Z  =  asymptotic  charge  of  the  atomic  core  (Z  s  Z  _  number  of  core 
electrons) , 

a  =  adjustable  constant. 

The  value  of  a  is  chosen  to  give  correct  asymptotic  behavior  of  P 

n  ^ 

This  form  of  V(r)  is  convenient  if  wave  functions  or  core  electrons  are  not 
known. 

III.  MHD  CODES 

Numerical  and  computational  support  and  conversion  to  VAX  is  provided 
for  the  cylindrical  core/corona  mode  "SPLAT"  and  the  Hartree-Fock  atomic 
physics  code  RCN  from  LASL  which  calculates  atomic  structure  parameters. 
The  SPLAT  code  is  a  time-dependent  self-similar  implosion  code  with  ID 
radial  fluid  grid  which  is  moved  in  a  Lagrangian  sense  and  enables  one  to 
track  the  local  neutral  plasma  position,  velocity,  and  acceleration  fields. 
The  ion  density  is  the  fundamental  fluid  quantity  which  is  evolved  and  the 
ion  fluid  is  heated  through  all  appropriate  physical  processes  and  cooled 
through  expansion  and  by  radiation.  The  electron  heating  takes  place  by 
compression,  viscous,  ohmic  heating  and  core/corona  heat  exchange  and 
radiative  losses.  The  RCN  code  from  LASL  is  used  for  the  self-consistent- 
field  calculation  of  atomic  radial  wave- functions  and  radial  integrals 
involved  in  the  calculation  of  atomic  energy  levels  and  spectra.  The 
calculation  is  of  radial  wavefunctions  of  a  spherically-symmetrized  atom, 
corresponding  to  the  center  of  gravity  energy  of  an  electron.  From  these 
radial  wavefunctions  the  coulomb  integrals  and  the  spin-orbit  integrals  are 
computed. 
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When  the  radial  component  F  .  of  the  eigenfunction  v  =  0_  .  R„  „ 

n,£  d  m,&  n,Jt 

is  multiplied  by  R#n  ^  then  the  product  Rr  gR*n  ^  gives  the  probability- 
density  distribution  T9*  as  a  function  of  the  nuclear-electron  distance  r. 


The  radial  equation 


13  (  2  3R  \ 
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^t-1  ■  R  +  (W-V)  R  =  0 


where  W  =  ± 


2ir2ueV 

n2h2 


V  = 


e2Z 


r  = 


nh 


&n2ye2Z 


has  a  solution 


n,JL 


H(n-H-I)!  Z3  \*  (2Zr\l  (  Zr_\  .  211+1  (zZr\ 
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and  the  expression  L  ntl^na^")  are  derivatives  of  the  Laquerre  polynomials. 


2  2 

^irr  (R  0)  is  the  probability-density  function  and 


The  function  D 

represents  the  probability  of  finding  the  electron  with  quantum  numbers  n,£ 
between  two  spheres  of  radius  r  and  r+dr.  Modeling  the  hydrogenic  electron 


we  can  plot  the  radial  function  R. 


(Rn,*> 


and  the  probability-density 


n,£  ’ 

function  D  as  a  function  of  the  electron  distance  r  for  various  atomic 
numbers  and  n  and 


Special  codes  and  models  were  developed  in  order  to  obtain  the  radial 
function  R  the  probability-density  factor,  and  the  density  distribution 
D  of  a  wave  function  of  an  electron.  Examples  of  these  graphs  are  given 
belcw  (Figures  1-9).  In  some  of  these  graphs  charge  distribution  of  non- 
hydrogenic  cases  are  compared  with  wave  functions  of  more  complicated 
programs . 

In  the  Bohr-Sommerfeld  calculations  of  the  elliptic  orbits  an  electron 
with  total  quantum  number  n  and  azimuthal  quantum  number  k  has  semimajor 
axis  a  and  semiminor  axis  b  given  by: 
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n  nk 

a:ai  r  *  b  =  “i  t 

where  is  the  first  Bohr  orbit  radius.  For  k  =  n  we  have  circular  orbits 
and  for  •  (n-l)  elliptic  orbits.  The  k=0  straight  line  orbit  is 

excluded.  The  discreteness  of  the  energy  states  in  the  Bohr-Sommerfeld 
case  is  introduced  quite  artificially.  On  the  contrary  in  the  case  of  the 
solution  of  the  Schrodinger  equation  the  discreteness  app  r  quite 
naturally  as  the  solution  of  the  Serodinger  equation. 

Tn  Figures  1  and  2  the  electron  charge  distribution  'on  of 

Cowan’s  LASL  calculations  for  A  2,  XI  2s  and  Is  electrons  is  co.  ,d  with 
our  (dotted  line)  calculation.  In  Figure  3  our  calculations  are  compared 
with  TRANP  (Milan)  calculations  for  the  A 2  XI  3d  electron.  The  units  of 
length  in  all  plots  are  in  first  Bohr  circular  orbit  radius. 

In  Figures  *1,  5,  and  6  the  Radial  factor  R  ,  the  probability  density 

p  2  30  p 

function  D  s  Jiur  (R  )  ,  and  the  Radial  factor  (R  „)  are  plotted 

^  IL  n 

respectively  in  3-D  for  the  hydrogenic  electron  n=8  and  J^3.  The  number  of 
peaks  (or  rings)  in  the  (R  )2  and  D  =  U-n r2  (R^)2  case  is  equal  to  the 
difference  (n-2).  For  the  n=8,  i?3  case  the  number  of  rings  is  five  and 
the  plots  show  the  radial  distribution  of  the  electron  around  the  nucleus 
with  five  maxima.  Since  we  consider  only  the  radial  component  Rn^  there  is 
a  rotation  symmetry  for  the  azimuthal  angle  6  and  the  angle  <J>. 

In  Figure  7  we  see  a  3-D  plot  of  the  probability-density  function  for 
the  hydrogenic  electron  2s  (n=2,  j?0).  In  the  cross  section  across  the  x 
axis  we  observe  two  maxima  for  the  peaks  since  n-2?2.  This  cross  section 
is  reminiscent  of  Figure  1.  There  we  had  the  charge  distribution  of  a  2s 
electron  in  A 2  XI.  There  the  peaks  were  at  different  distances  from  the 
center  of  the  nucleus.  Figure  8  shows  the  whole  view  of  Figure  7  from  a 
different  poslton  of  the  camera. 
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The  solution  of  the  differential  equation 
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+  £  (£,  +  1  )  =  0 


is  the  angular  component  ©  of  the  eigenfunction  y 

m  fJC 

/  (2£  +  1 )  (£-m) !  \*  m  pm  .  . 

0m,£  =  \  2(£+m) !  )  *  sin  e  *  P£  (e0Se) 


where  P®  (cose)  is  the  associated  Legendre  polynomial.  The  product  0m  ^ 
0*m  •  sinede  is  the  probability  of  finding  the  electron  between  e  and 

e+de- 


The  complete  eigenfunction  is 

M'_  „  _  =  4>  •  0_  „•  _ 

m,£  ,n  m  m,£  £  ,n 

and  the  probability-density  factors  (for  symmetric  <j>  )  for-  the  radial  and 

2  2 

angular  components  (©_  )  and  (R  )  give  the  probability  densities  for 

mj,  £  .n 

finding  the  electron  between  two  spheres  of  radius  r  and  r+dr  and  azimuthal 
angle  between  0  and  0+dQ. 

Programs  calculating  the  Legendre  polynomials  for  the  angular 
component  of  the  wave  function  are  developed.  Using  the  radial  and  angular 
distribution  the  probability  density  of  hydrogenic  electrons  of  any 
combination  of  n,  £,  and  m  can  be  plotted.  A  few  examples  are  given  in 
Figures  9  through  15. 

The  number  N  of  the  wings  in  the  butterfly  of  the  angular  component 
for  an  electron  with  quantum  numbers  £  and  m  is  given  by: 

N  :  2*  (£-m+1 ) 

In  Figures  10,  11,  12,  and  13  we  see  how  the  angular  distributions 
®33’  e32’  ©31’  and  ©30  chan8e  the  probability  density  of  the  hydrogenic 
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electron  n=8,  £.=3  in  the  five  rings  of  the  radial  distribution  Rg^  of 
Figure  9.  Figures  13,  I1)  and  15  show  the  radial  and  angular  probability- 
density  for  different  combinations  of  the  quantum  numbers  n,Sl,  and  m. 

The  old  SERDY  program  was  restructured  and  improved.  SERDY  is  a 
hydro-transport,  radiation  energy  deposition,  radiation  transport  and 
thermal  conduction  program.  The  physics  in  this  spherical-equilibrium- 
radiation-dynamics  program  with  respect  to  the  hydrodynamic  equations 
accounting  for  radiation  energy  pressure  and  radiant  heat  exchange  was 
corrected  so  that  an  investigation  of  an  argon  seeded  microballoon  implosion 
could  be  possible.  The  results  of  this  investigation  were  presented  at  the 
APS,  Division  of  Plasma  Physics  Meeting,  October  1981.  The  following  are 
abstracts  of  two  papers  delivered  at  the  APS,  New  York  Meeting. 

CRE  Analysis  of  Argon  Puff  Gas  R-Series  Emission  Spectra*. 

K.G.  WHITNEY,  C.  AGRITELLIS*,  AND  J.  DAVIS,  Naval  Research  Laboratory.  —In 
collisional-radiative  equilibrium  (CRE),  the  strength  of  argon  puff  gas 
spectral  features,  both  relative  and  absolute,  is  a  direct  reflection  of 
the  distribution  of  strengths  of  collisional  and  radiative  couplings  within 
the  plasma.  A  detailed  atomic  model  of  argon  has  been  constructed 
containing  extensive  level  structure  in  Ar  XVI,  XVII,  and  XVIII,  Lyman  and 
Balmer  series  line  couplings  in  Ar  XVII  and  Ar  XVIII,  eight  of  the  dominant 
line  couplings  in  Ar  XVI,  and  Gaunt  factor  contributions  to  the  continuum 
free-bound  emission  background.  Using  this  model,  we  will  describe  the 
behavior  of  instrumentally  broadened,  optically  thick  argon  K-shell 
emission  spectra  under  the  influence  of  temperature  and  density  gradients 
and  as  a  function  of  the  degree  of  plasma  compression  and  plasma  mass. 

+  This  work  was  supported  by  DNA 
*  Science  Applications,  Inc.  McLean,  Va. 

Radiation  Hydrodynamics  of  a  Radiation  Imploded  Argon  Microballoon*. 

C.  AGRITELLIS*,  K.  G.  WHITNEY  and  J.  DAVIS,  Naval  Research  Laboratory.  — 
The  radiation  hydrodynamics  of  a  radiation  imploded  argon  microballoon  is 
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investigated.  The  mode  includes  an  absorption  mechanism,  non-LTE  radiation 
physics,  line  and  continuum  radiation  transfer,  spectral  line  profile 
functions  and  hydrodynamics.  The  purpose  of  this  investigation  is  to  study 

the  behavior  of  detailed  K-shell  emission  spectra  as  a  function  of  initial 
driver  conditions.  Detailed  comparisons  will  be  made  between  the  continuum 
intensity  spectral  line  intensity  and  widths,  and  continuum  edge  effects 
for  several  initial  conditions. 

+  Work  supported  by  Defense  Nuclear  Agency 
*  Science  Applications,  Inc.  McLean,  Va. 

While  X-ray  spectra  from  imploded  plasmas  are  often  used  to  infer  an 
average  plasma  temperature  and  density,  these  are  necessarily  uneconomical 
diagnostics  since  there  is  generally  a  considerable  amount  of  space-time 
variability  to  plasma  temperatures  and  densities  during  the  compression  and 
subsequent  expansion  stages.  We  have  investigated  the  influence  of 
temperature  and  density  gradients,  at  an  instant  of  time,  on  the  structure 
of  the  collisional-radiative  couplings  within  a  compressing,  spherical 
plasma  and  correlated  these  structures  with  the  shape  and  intensity  of  the 
Instantaneous  emission  spectrum.  It  is  theoretically  possible  to 
accurately  generate  and  analyze  temporally  and  spectrally  resolved 
broadband  emission  spectra  from  imploded  plasmas  that  leave  complicated 
temperature  and  density  profiles.  It  is  important  to  obtain  absolute 
(spectrally  resolved)  calorimeter  data,  plasma  size  information,  as  well  as 
film  spectral  data  in  order  to  make  unambiguous  theoretical  interpretations 
of  emission  spectra  from  optically  thick  plasmas.  The  so-called  small 
discrepancies  that  are  seen  in  line  ratios  that  are  used  to  predict  the 
temperature  for  imploded  plasmas  may  be  indicators  that  the  plasma  has  no 
well  defined  temperature  but  rather  large  temperature  and  density 
gradients.  For  each  of  the  argon  emission  features  within  the  spectrum,  we 
use  the  ratios  of  photo-excitation  to  electron  collisional  excitation  as  a 

4 

measure  of  the  structure  of  the  collisional-radiative  couplings  within  the 
plasma.  These  ratios,  which  vary  throughout  the  plasma  volume,  are  plotted 
in  Figure  17  and  18  in  the  form  of  a  histogram  with  plasma  cells 
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represented  along  the  y-axis  and  the  lines  and  edges  represented  along  the 
|  x-axis.  Low  frequency  lines  are  characterized  by  high  electron  collisional 

excitation  values  while  the  high  frequency  lines  and  the  edges  have  very 
high  photo  excitation  performance. 

|  Graphic  programs  and  packages  are  created  for  diagnostic  purposes  and 

multi-color  presentation  of  non-LTE  radiation  from  argon  seeded 
microballoon  implosion  using  NRL/LASNEX  and  NRL  hydrodynamics  models. 
Color  pictures  were  created  corresponding  to  K-Shell  X-Ray  emission  in  the 
|  compression  phase  from  an  argon  microballoon  using  the  NRL/LASNEX  radiation 

model.  In  Figures  19  and  20  the  electron  temperature  and  density  gradients 
are  shown  in  a  compressed  spherical  plasma  along  with  the  shape  and 
intensity  of  the  instantaneous  emission  spectrum.  TVo  color  movies  of  the 
|  compression  and  a  set  of  pictures  of  the  full  spectrum  at  different  times 

provide  a  spectacular  view  of  the  whole  event.  A  color  movie  for  the  NRL 
radiation  case  and  a  set  of  color  pictures  for  both  APS  presented  papers 
have  already  been  delivered  to  NRL. 
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Fig.  1.  Charge  Distribution  of  a  2s  electron  in  Al  XI  Compared  with 
LASL  (Covan)  Program 
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Fig.  2.  Charge  Distribution  of  a  Is  electron  in  Ai  XI  Compared  with 
LASL  (Cowan)  Program 
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CHRRGE  DISTRIBUTION  FOR  RL  XI 


Fig.  3.  Charge  Distribution  of  a  3d  electron  in  Ai  XI  Compared  with 
NRL-MILAN  Program 
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Fig.  4.  Radial  Factor  R^  £  for  Hydrogen  n  -  8,  l  -  3 
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Fig.  8.  Probability-Density  Function  D 
2s  Electron 
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Fig.  19.  Black  and  White  Copy  from  the  Color  Movie  of  the  NRL/LASNEX 
Case 
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